Chaplygin inflation in loop quantum cosmology 
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In this paper we discuss the inflationary universe in the context of a Chaplygin gas equation of 
state within the framework of the effective theory of loop quantum cosmology. Under the slow- 
roll approximation, we calculate the primordial perturbations for this model. We give the general 
expressions of the scalar spectral index, its running, and the tensor-to-scalar ratio, etc. For the 
chaotic inflation with a quadratic potential, using the WMAP 5-year results, we determine the 
parameters of the Chaplygin inflation model in loop quantum cosmology. The results are consistent 
with the WMAP observations. 



I. INTRODUCTION 

The inflation paradigm provides a compelling solu- 
tion for many long-standing problems of the cosmologi- 
cal standard model (such as flatness, horizon, monopoles, 
etc.) by positing an epoch of accelerated expansion in the 
early universe [1-5]. This accelerated period of expan- 
sion also generates superhorizon fluctuations and thus 
predicts an almost scale-invariant density perturbation 
power spectrum, which has received strong observational 
support from the measurement of the temperature fluc- 
tuation in the cosmic microwave background (CMB) ra- 
diation [6-12]. Conceptually, however, the inflationary 
scenario is incomplete due to the existence of the big 
bang singularity [13]. Einstein's classical theory of gen- 
eral relativity (GR) breaks down near such a singularity 
since quantum effects are expected to be important at 
very high energies in the early universe. Thus, the classi- 
cal theory of GR has to be replaced by some theoretical 
framework of quantum gravity which should remain well 
defined even at very high curvatures. 

Loop quantum gravity (LQG) is a leading nonpertur- 
bative background independent approach to quantizing 
gravity [14-16]. The underlying geometry in LQG is 
discrete and the continuum spacetime can be obtained 
from the quantum geometry in a large eigenvalue limit. 
Loop quantum cosmology (LQG) focuses on symmetry 
reduced models (with homogeneous and isotropic space) 
but inherits quantization scheme and techniques from 
LQG [17]. Within the framework of LQC, some long- 
standing issues concerning the quantum nature of the 
big-bang are resolved in the context of homogeneous and 
isotropic universe with a scalar field. Using extensive an- 
alytical and numerical methods, the analysis of the evolu- 
tion of the semiclassical states for a spatially flat universe 
has shown that the universe has a pre-big-bang branch, 
joined deterministically to the post-big-bang branch by a 
quantum bounce in the deep Planck regime through the 
LQG evolution [18, 20, 21]. Thanks to the nonperturba- 
tive background independent methods of LQG, the idea 
of the nonsingular bounce can be realized in a natural 
fashion. 

Recently, it has been shown that the discrete quantum 



dynamics can be well approximated by an effective mod- 
ified Friedmann dynamics [20-22]. The modified Fried- 
mann equation can be obtained from the effective Hamil- 
tonian constraint, which can be used to investigate the 
role of nonperturbative quantum correction conveniently. 
It is remarkable that the quantum geometric effects lead 
to a modification to the Friedmann equation at the 
scales when p becomes comparable to a critical density 
Pc which is close to the Planck density {pc ~ 0.826*^^) 
[20, 21, 23-25]. Within the framework of LQC, the infla- 
tionary universe model has been studied in detail [26] . 

On the other hand, it is well known that the general- 
ized Chaplygin gas (GGG) is an alternative candidate for 
explaining the acceleration of the universe. The GGG is 
described by an exotic equation of state in the form [27] 
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where pch and pch are, respectively, the energy density 
and pressure of the GCG, a is a constant in the range 
< a ^ 1, and A is a positive constant. It should also 
be mentioned, however, that "superluminal" values a > 1 
are well possible too, as was recently shown in Ref. [28]; 
causality requires only some change in Eq. (1) for pch 
very close to A^^". The case a — I corresponds to the 
original Chaplygin gas [29, 30]. Inserting this equation of 
state into the stress-energy conservation equation leads 
to the following energy density 
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where a is the scale factor of the universe and _B is a 
positive integration constant. The model of GCG as well 
as its further generalization have been extensively studied 
in the literature [31-42]. In addition, it should be pointed 
out that the GCG model may be viewed as a modification 
of gravity [43] . 

Recently, it was considered that the GCG model can 
also be used to describe the early universe. In this con- 
sideration, Eq. (2) is not viewed as a consequence of 
the equation of state (1), but rather, as arising due to a 
modification of gravity itself. Furthermore, one can also 
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assume that during inflation the gravity dynamics may 
give rise to a modified Friedmann equation [44] 



A + 



(3) 



where Mp is the reduced Planck mass 



is the energy 

density of the inflaton field. This is the so-called Chaply- 
gin inspired inflation scenario [44] , since the modification 
in Eq. (3) is realized firom an extrapolation of Eq. (2): 
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(4) 



corresponds to the matter energy den- 



where Pm ^ 

sity. Following this idea, lately, some work has been done, 
this involves tachyon-Chaplygin inflation [45], warm- 
Chaplygin inflation [46] , Chaplygin inflation on a brane 
[47], and so forth [48]. In this paper, we shall study the 
Chaplygin inspired inflationary model in loop quantum 
cosmology. 

This paper is organized as follows. In the next section 
we briefly review the efl'ective theory of LQC. In section 
II, we describe the Chaplygin inflationary universe in the 
framework of LQC. In section IV, we study the primor- 
dial perturbations and give the expressions of scalar and 
tensor perturbations for this model. Also, we shall com- 
pare our results to the WMAP five-year observations. 
Conclusion is given in section V. 



II. EFFECTIVE DYNAMICS IN LOOP 
QUANTUM COSMOLOGY 

LQG is a canonical quantization of gravity based upon 
Ashtekar-Barbero connection variables. The phase space 
of classical GR in LQG is spanned by SU(2) connection 
A\ and the triad E'f on a 3-manifold M. (labels a and 
i denote space and internal indices respectively), which 
are two conjugate variables encoding curvature and spa- 
tial geometry, respectively. Likewise, LQC is a canoni- 
cal quantization of homogeneous spacetimes based upon 
techniques used in LQG. In LQC, due to the symmetries 
of the homogeneous and isotropic spacetime, the phase 
space structure is simplified, i.e., the connection is deter- 
mined by a single quantity labeled c and likewise the triad 
is determined by a parameter p . The variables c and p are 
canonically conjugate with Poisson bracket {c, p} = 7/3, 
here we have used the unit Mp = 1 for convenience, and 7 
is the dimensionless Barbcro-Immirzi parameter which is 
set to be 7 w 0.2375 by the black hole thermodynamics in 
LQG. For the spatially flat model of cosmology, the new 
variables have the relations with the metric components 
of the Friedmann-Robertson- Walker (FRW) universe as 



c = 7a, p = a 



(5) 



where a is the scale factor of the universe. Classically in 
terms of the connection-triads variables the Hamiltonian 



constraint is given by 



-h ZtM, 



(6) 



where Hm is the matter Hamiltonian. 

The elementary variables used for quantization in LQC 

are the triads and holonomies of the connection. The 
holonomy over an edge of a loop is defined as = 
cos(/zc/2)-|-2sin(/ic/2)Ti, where is related to Pauli spin 
matrices as = — i(T,/2 and dimensionless p is related to 
the physical length of the edge over which holonomy is 
evaluated (note that p is also the eigenvalue of the triad 
operator p). In the Hamiltonian formulation for homoge- 
neous and isotropic spacetime, the dynamical equations 
can be determined completely by the Hamiltonian con- 
straint. Under quantization, the Hamiltonian constraint 
gets promoted to an operator and the quantum wave 
functions are annihilated by the operator of the Hamil- 
tonian constraint. In LQC, it is expected that modifica- 
tions due to LQC effects will appear in the Hamiltonian 
constraint, and from the modified Hamiltonian constraint 
the effective Friedmann constraint will be derived. In 
quantization the Hamiltonian constraint operator is ob- 
tained by promoting the holonomies and the triads to 
the corresponding operators. Consequently, this leads to 
a discrete quantum difference equation, which indicates 
that the underlying geometry in LQC is discrete [19]. In- 
terestingly, the solutions of this difference equation are 
nonsingular. 

So far we see that the underlying dynamics in LQC 
is governed by a discrete quantum difference equation in 
quantum geometry. However, an effective Hamiltonian 
description on a continuum spacetime can be constructed 
by using semiclassical states, which has been shown to 
very well approximate the quantum dynamics [20, 21]. 
This analysis reveals that on backward evolution of our 
expanding phase of the universe, the universe bounces at 
a critical density (near the big bang singularity) into a 
contracting branch [18, 25]. Thus the classical singular 
problem can be successfully overcome within the context 
of LQC by a nonsingular bounce. In addition, the effec- 
tive equations for the modified Friedmann dynamics can 
be derived from the effective Hamiltonian constraint with 
loop quantum modifications, which can be used to inves- 
tigate the role of nonperturbative quantum corrections. 
An important feature for the modified dynamics is that 
a p^ term which is relevant in the high energy regime is 
included in the classical Friedmann equation. The modi- 
fied term is negative definite implying a bounce when the 
energy density reaches a critical value on the order of the 
Planck density. 

The effective Hamiltonian constraint, to leading order, 
is given by [22] 



n 



eff 



asm^{pc) + Hm, 



(7) 



where p is the kinematical length of the edge of a square 
loop which has the area given by the minimum eigenvalue 
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of the area operator in LQG; the area is ^ = /i^a^ = al"^, 
where a is of the order unity and Ip is the Planck length. 
The modified Friedmann equation can then be derived 
by using the Hamilton's equation for p, 

p = {p,Weff} = -g^i— = — sm(;Uc)cos(Mc), (8) 

which combined with Eq. (5) yields the rate of change of 
the scale factor 



a= sin(/xc) cos(^c). 
7/x 



(9) 



Furthermore, the vanishing of the Hamiltonian con- 
straint, Hen ~ 0, implies 



sin^(/xc) 



3a 



(10) 



Combining Eqs. (9) and (10) yields the effective Fried- 
mann equation for the Hubble rate H = a/ a, 



1 



1 - — 
3^V^ Pc 

with the critical density given by 
Pc = 4737-3. 



(11) 



(12) 



The modified Friedmann equation provides an effective 
description for LQC which very well approximates the 
underlying discrete quantum dynamics. The nonpertur- 
bative quantum geometric effects are manifested in the 
modified Friedmann equation with a correction term. 
The negative definition of the term implies that the 
Hubble parameter vanishes when p = p^ and the uni- 
verse experiences a t\irnaro\.md in the scale factor. When 
p Pc, the modifications to the Friedmann equation be- 
come negligible, and the standard Friedmann equation 
is recovered. In addition, it should be noted that, inter- 
estingly, p^ modifications also appear in string inspired 
braneworld scenarios and it has been shown that there ex- 
ist interesting dualities between the two frameworks [23] . 
Such modifications in braneworlds, however, are usually 
positive definite so that a bounce is absent, unless the ex- 
istence of two timelike extra dimensions is assumed [49] . 
For extensive studies in this framework, see e.g. [50-57]. 
Other mechanism of realizing a bouncing (or oscillating) 
universe can be found in e.g. [58-63]. 



III. CHAPLYGIN INFLATIONARY UNIVERSE 
IN LOOP QUANTUM COSMOLOGY 

In this section, we shall investigate Chaplygin inspired 
inflation within the framework of the effective theory of 
LQC. In this case the dynamics of the inflationary uni- 
verse is governed by the following equation. 



3F^=(A + p(^+"))-° 



1 - 



. (13) 



Here, p^ is the canonical inflaton, p^ = (j)/2 + V{4>), and 
V{(j)) is its potential. It should be noted that in the 
low-energy regime, [A + p^^"'""^]i/(i+«) <^ p^, the stan- 
dard Chaplygin inflation model is recovered; while in the 
high-energy regime, the quantum gravity effects leads to 
a super-inflation that is in the "primary inflation" . How- 
ever, what is of interest for us is the "observable inflation" 
with the last 60 e- foldings. We thus only discuss the sub- 
sequent normal inflation stage (after the super-inflation 
phase) in what follows. Also, in the following we will take 
a = 1 for simplicity, which means the usual Chaplygin 
gas. 

For the inflaton fleld, its equation of motion is of the 
standard form. 



3ff0 + V' = Q, 



(14) 



where dots denote the derivatives with respect to the cos- 
mic time and V = dV{(l))/d(l), and the Hubble parameter 
H is given by Eq. (13). If 0^ < v{(j)) and 4> < 3H<j), 
the scalar fleld will slowly roll down its potential, and 
the exact evolution equation (14) can be replaced by the 
slow-roll approximation 



-V'/3H. 



(15) 



Under the slow-roll approximation, the energy density of 
the scalar field approximates as p ~ y{4>)y the Friedmann 
equation (13) consequently can be written as 



We can also define the slow-roll parameters 



H 



1 VV'2 (1 - 



2^/A+V^^ 



2(^ + 1/2)3/2 



Pc > 



77 = e + 5 : 



V" 



(A + y2)i/2 



(16) 



(17) 



(18) 



where 5 = —(j)/H<j), by definition. The slow-roll condition 
can be expressed as e, jvyj <C 1. 

Note that in the limit A ^ 0, the slow-roll parame- 
ters e and rj coincide with the case in inflation model in 

LQC [26]. Also, in the low-energy limit, 



/>1 + p; < Pc, 

the slow-roll parameters reduce to the standard form of 
Chaplygin inflation. 

The parameter e quantifles how much the Hubble rate 
H changes with time during inflation. Inflation can be 

attained only if e < 1. This condition could be writ- 
ten explicitly in terms of the scalar potential V and its 
derivative V 



Pc 



VV'1 

< (^ + ^2)3/2 



Pc 



(19) 
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Inflation ends when e < 1 ceases to be satisfied, which 
implies 



VfV}2 



1 



= (A + y/)3/2 



p. 



2-.1/2 



(20) 



where the subscript / is used to mark the end of inflation. 

The small change of e-foldings satisfies dN = —Hdt; 
using Eqs. (15) and (16), one obtains the number of e- 
foldings of slow-roll inflation remaining at a given epoch, 



^(<^*) = J 

<t>j 



V' 



1- 



d(j), 



(21) 



where the subscript * marks the epoch when the cosmo- 
logical scales exit the horizon. 



IV. PRIMORDIAL PERTURBATIONS AND 
WMAP FIVE-YEAR RESULTS 

In this section we shall study the primordial perturba- 
tions of this inflationary model. We shall give the expres- 
sions of scalar and tensor perturbations for the Chaplygin 
inflation in LQC. 

The perturbation 5(j) can be treated as a massless free 
field, and its vacuum fluctuation can be regarded as a 
classical quantity a few Hubble times after horizon exit. 
The spectrum of the field perturbation is 



(22) 



The corresponding curvature perturbation is given by 
n = {H/^)5(l). Using Eqs. (15), (16) and (22), we get 
the amplitude of scalar perturbation as 



1 - 



(23) 

which is evaluated at the Hubble radius crossing k = aH. 
Note that in the limit ^ ^ 0, the amplitude of the scalar 
perturbation given by Eq. (23) coincides with Rcf. [26]. 

Since H is slowly varying, we have the relation 
d\nk{(j)) = -dN{4>). Then, from Eq. (23) we get the 
scalar spectral index, n^ — 1 = dlnA^/rflnfc = —&e + 2ri, 
or equivalently, 




(24) 



Once again, we notice that in the limit A^ Q the scalar 
spectral index coincides with that in Ref. [26]. 



Another important quantity for inflationary model is 
the running of the scalar spectral index = diis/dlnk. 
Prom Eq. (24), we can get the running of the scalar 
spectral index for this model 



AjA + V^) 



1 - 



Pc 



1 - 



[3e,^-??,^]e, (25) 



where e and rj are given by Eqs. (17) and (18), and 
the subscript ,0 labels the derivative with respect to the 
inflaton field <j>. 

Inflation also generates gravitational waves with two 
independent components /i+,x which have the same ac- 
tion as a massless scalar field. Likewise, the amplitude 
of tensor perturbation can also be given. 



27r J 



37r2 I Pc 



(26) 

which is also evaluated at the horizon exit k = aH. 
The spectral index of tensor perturbation is given by 
rit = dAf^/dlnk that will not be written explicitly here. 
Then, from the expressions (23) and (26), we can write 
the tensor-to-scalar ratio, at fc = aH, as 



8V' 



12 



{A + v^)[i-{A + v^y/ypc 



(27) 



Recently, the Wilkinson Microwave Anisotropy Probe 
(WMAP) five-year data were released [12]. The WMAP 
5-year data provide stringent constraints on inflationary 
models. If we assume that the primordial fluctuations 
are adiabatic with a power law spectrum, the WMAP 
5-year data (WMAP only) give: = 0.9631°;^^^; while 
the WMAP data combined with the distance measure- 
ments from the Type la supernovae (SN) and Baryon 
Acoustic Oscillations (BAO) in the distribution of galax- 
ies (WMAP-KBAO-FSN) give: = 0.960tggJt The 
amplitude of curvature perturbation is = 2.4 X 10"^, 
measured by WMAP at fc* = 0.002 Mpc"^ With the 
WMAP data combined with BAO and SN, it is found 
that the limit on the tensor-to-scalar ratio is r < 0.20 
(95% CL) and rig > 1 is disfavored even when gravi- 
tational waves are included. We will make use of the 
WMAP 5-year results to determine the parameters of 
the model. 

Consider now an inflaton scalar field 4> with a chaotic 
potential. We focus on a quadratic potential V{(j)) = 
m^0^/2, where m is the mass of the scalar field. First, 
we should evaluate (f>^ (that corresponds to the time of 
horizon-crossing) by using Eq. (21). In view of that the 
LQC effect is not prominent (for convenience we define 
1/ = (A + V'^) / Pc and we have u < 10~^, see [26]), we ne- 
glect 1/ in this evaluation. Then, the integral in Eq. (21) 
can be given directly (see Eq. (19) in [44]). Furthermore, 
considering 0/^0 and the unimportance of logarithm 
term, one obtains </)« = 2-\/iV [44]. 
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Next, using the WMAP 5-year results 
(WMAP+BAO+SN), A^j = 2.4 x 10"^ and = 0.960, 
we can determine the parameters of the LQC Chap- 
lygin inflation model. Taking N = 55, we derive 
A = 3.36 X 10-18 and m = 5.9 x 10-^ note that 
here the unit Mp = 1 has been used. Based on these 
values of parameters, we can further evaluate the other 
observational quantities. For example, we give the 
nmning of the scalar spectral index as = —6.8 x 10"^ 
and the tensor-to-scalar ratio r = 0.1, which are in good 
consistence with the WMAP observations. Also, we 
verified the fact that the LQC effect can only leads to 
a very tiny imprint in the primordial power spectrum, 
1/ = 8.5 X 10-1°, which is in good accordance with the 
result in Ref. [26]. 

V. CONCLUSION 

In this paper we have studied the Chaplygin inflation 
model in the framework of the effective theory of LQC. 
In LQC, the nonperturbative quantum geometry effects 
lead to a term with a negative sign in the modified 
Friedmann equation. On the other hand, inspired by the 
Chaplygin equation of state, a phenomenological mod- 
ification of gravity is also considered in the Friedmann 
dynamics. The LQC Chaplygin inflationary model com- 
bines the above both features. 

Under the slow-roll approximation, we calculated the 
primordial perturbations for this model. We gave the 
general expressions of the scalar spectral index Ug, its 
running as, and the tensor-to-scalar ratio r, etc. For 



the chaotic inflation with a quadratic potential, using 
the WMAP 5-year results, A|,(A:») = 2.4 x 10"^ and 
ns{k.^,) = 0.960, we determined the parameters of the 
LQC Chaplygin inflation model, A = 3.36 x IQ-i* and 
m = 5.9 X 10-^. This leads to the results of the running of 
the scalar spectral index as(A;*) = —6.8 x IQ-^ and the 
tensor-to-scalar ratio r(fc*) = 0.1, which are consistent 
with the WMAP observations. 

The theory of LQC gives rise to a quantum bounce in 
the high energy regime when the loop quantum effects are 
dominative power. The classical big-bang singularity is 
thus replaced by the quantum bounce. After the bounce, 
a super-inflation phase was emergent in a natural way. 
Then the universe underwent a normal inflation stage. 
In this paper, we restrict that the Chaplygin-inflation 
happens in the normal inflation phase. This leads to that 
the imprint of the LQC effect in CMB sky is too weak 
to be observed by present observational data (the LQC 
parameter ;/ < 10-^, sec Ref. [26]). For the Chaplygin 
inflation case with quadratic potential, we derived v = 
8.5 X 10-1° that is in agreement with the conclusion of 
Ref. [26]. If we set the Chaplygin inflation to happen 
in the super-inflation phase, it is believed that the LQC 
effect should be prominent and the imprint of LQC effect 
may be detected in the CMB sky. We hope to return to 
this point in the near future. 
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